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ABSTRACT 

We present deep near-infrared JHK imaging of four 10' x 10' fields. The observations were carried 
out as part of the Multiwavelength Survey by Yale-Chile (MUSYC) with ISPI on the CTIO 4m 
telescope. The typical point source limiting depths are J ~ 22.5, H ~ 21.5, and K ~ 21 (5a; Vega). 
The effective seeing in the final images is ~ l'.'O. We combine these data with MUSYC UBVRIz 
imaging to create iGselected catalogs that are unique for their uniform size, depth, filter coverage, 
and image quality. We investigate the rest-frame optical colors and photometric redshifts of galaxies 
that are selected using common color selection techniques, including distant red galaxies (DRGs), 
star- forming and passive BzKs, and the rest-frame UV-selected BM, BX, and Lyman break galaxies 
(LBGs). These techniques are effective at isolating large samples of high redshift galaxies, but none 
provide complete or uniform samples across the targeted redshift ranges. The DRG and BM/BX/LBG 
criteria identify populations of red and blue galaxies, respectively, as they were designed to do. The 
star-forming BzKs have a very wide redshift distribution, a wide range of colors, and may include 
galaxies with very low specific star formation rates. In comparison, the passive BzKs are fewer 
in number, have a different distribution of K magnitudes, and have a somewhat different redshift 
distribution. By combining these color selection criteria, it appears possible to define a reasonably 
complete sample of galaxies to our flux limit over specific redshift ranges. However, the redshift 
dependence of both the completeness and sampled range of rest-frame colors poses an ultimate limit 
to the usefulness of these techniques. 

Subject headings: catalogs — surveys — galaxies: distances and redshifts — galaxies: high-redshift 
— infrared: galaxies 



1. INTRODUCTION 

The issues surrounding galaxy formation and evolution 
provide significant motivation for ongoing astrophysical 
research. Although several of the basic physical processes 
that drive the formation and evolution of galaxies were 
identified some time ago, such as the gravitational col- 
lapse of dark matter, the cooling and dissipative collapse 
of baryons, the variation in star formation rates with 
time, and galaxy mergers, constructing detailed models 
that reproduce observed properties for a wide range of 
galaxy types and redshifts has proven exceedingly dif- 
ficult. A precise determination of the cosmological pa- 
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rameters has removed some of the uncertainties, but it 
is largely accepted that comprehensive observations of 
galaxies at low and high redshifts are necessary to refine 
our understanding of galaxy evolution. 

Recent years have seen dramatic progress in the obser- 
vational study of high redshift galaxies. Much of this 
progress has been driven by deep imaging of "blank" 
fields with multiple bandpasses. These surveys often rely 
on carefully-designed color selection techniques, mak- 
ing use of only a few bands, to isolate samples of high 
redshift galaxies. The most well known color selection 
criteria identify the Lyman Break Galaxies (LBGs) at 
z ~ 3 (Steidel et al. 1996), using the U n GR bands. This 
technique has been extended to select LBGs at higher 
redshifts, and to select similar star-forming galaxies at 
z ~ 2.3 and z ~ 1.7 (BX and BM galaxies, respectively; 
Adelberger et al. 2004). Selecting galaxies in optical 
bands may miss many red high redshift galaxies, so NIR 
selection techniques have also become important. Ex- 
tremely Red Objects (EROs) are at z > 1, and have been 
selected using a variety of criteria, such as R — K > 5 
and / - H > 3 (McCarthy 2004). Franx et al. (2003) 
used J — K > 2.3 to select distant red galaxies (DRGs) 
at z ~ 2 — 4. Daddi et al. (2004) proposed a technique 
involving the BzK bands to identify galaxies at z > 1.4. 
Yan et al. (2004) used (z — 3.6^iti)ab > 3.25 to isolate 
IRAC Extremely Red Objects (IEROs). 

The primary advantage of such selection techniques is 
that they rely on only a few observed bands to isolate 
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large samples of galaxies in a (hopefully) well-defined red- 
shift range. However, diagnostic information about indi- 
vidual galaxies or about the range of properties spanned 
by a sample of galaxies requires more bands to trace the 
detailed spectral energy distributions over a wide wave- 
length range. Specific color selection techniques may also 
identify only a subset of galaxies in the targeted redshift 
range (by design), offering a limited view of the diver- 
sity of galaxy properties at that redshift. Surveys with 
a large number of observed bands may be able to use 
photometric redshifts to obtain a more complete sample 
of galaxies in any given redshift range, and are better 
able to specify the intrinsic properties of high redshift 
galaxies. 

With these issues in mind, we have executed a deep op- 
tical/NIR survey of four southern and equatorial fields 
as part of the Multiwavclcngth Survey by Yale-Chile 
(MUSYC). 13 Subsets of these data have been used to 
study several characteristics of if-sclcctcd galaxies at 
z > 2, including the number density and colors of mas- 
sive galaxies (van Dokkum et al. 2006), the clustering 
properties (Quadri et al. 2006), and the luminosity func- 
tion (Marchesini et al. 2006). Additionally, Kriek et al. 
(2006a, b,c) obtained NIR spectroscopy of the brightest 
MUSYC galaxies, and Webb et al. (2006) used Spitzer 
MIPS imaging to infer the star formation rate of MUSYC 
DRGs. 

The purpose of this paper is two-fold. First, we present 
the MUSYC deep NIR imaging. We describe the ob- 
servations, data reduction, and the if-selected catalogs 
with full UBVRIzJHK photometry. Second, we use 
these catalogs and high quality photometric redshifts to 
investigate the properties of galaxies that are selected 
with common color selection criteria. We use Vega mag- 
nitudes unless noted otherwise. Throughout we use 
H = 70 km^s^Mpc- 1 , Q m = 0.3, and ft A = 0.7. 

2. THE MULTIWAVELENGTH SURVEY BY YALE-CHILE 

The MUSYC survey consists of several components: 
optical imaging of four 30' x 30' fields, NIR imag- 
ing over the same area, deeper NIR imaging over four 
10' x 10' subfields, and spectroscopic followup. Some of 
the MUSYC fields also benefit from imaging by Chandra 
and XMM in the x-ray, HST in the optical, and Spitzer 
in the infrared. 

The survey design is described by Gawiser et al. 
(2006a). The four 30' x 30' fields were selected to have 
low Galactic reddening, HI column density (Burstein & 
Heiles 1978), and 100/im dust emission (Schlegel et al. 
1998). They were also chosen to have high Galactic lat- 
itude to reduce the number of stars, to cover a wide 
range in right ascension to enable flexible scheduling, 
and to be accessible from Chilean observatories. These 
fields were observed with UBVRIz filters for complete 
optical coverage. A narrowband 5000A filter was also 
used to identify Lyman a emitters at z ~ 3 (Gawiser et 
al. 2006b). Gawiser et al. (2006a) presents images and 
optically-selected catalogs of one of the MUSYC fields, 
EHDF-S. The full optical data for the remaining fields 
will be described elsewhere (E. Gawiser et al. 2007, in 
preparation) . 

The four 30' x 30' fields were also observed in at least 

13 http://www.astro.yale.edu/MUSYC 



one of the J, H, and K bands. We refer to this as the 
"wide" portion of the MUSYC NIR imaging. These data 
were collected and processed in a similar way to the data 
described in this paper (E. Taylor et al. 2007, in prepa- 
ration; G. Blanc 2007 et al. 2007, in preparation). The 
typical depths are J ~ 22. and K ~ 20. 

Four 10' x 10' subfields were observed to greater depth 
in all of JHK. These data, which are referred to as the 
"deep" NIR MUSYC imaging, are presented in this pa- 
per. Two of these subfields, HDFS1 and HDFS2, are 
adjacent and lie within the larger 30' x 30' MUSYC Ex- 
tended Hubble Deep Field-South (EHDF-S). The other 
two subfields lie within the larger MUSYC 1030 and 1255 
fields. We did not perform deep NIR imaging in the 
fourth large MUSYC field, ECDF-S, because very deep 
imaging in the central region of this field has been made 
available by the Great Observatories Origins Deep Sur- 
vey (GOODS) team (Giavalisco et al. 2004). The exact 
locations of the deep MUSYC fields within the larger 
30' x 30' fields were chosen to avoid bright stars. Infor- 
mation about these fields is given in Table 1. 

3. OBSERVATIONS 

The MUSYC deep NIR observations were performed 
with the Infrared Sideport Imager (ISPI) (Probst et al. 
2003; van der Bliek et al. 2004) on the CTIO Blanco 
4m Telescope. The detector is a 2048 x 2048 HgCdTc 
HAWAII-2 array, with a ~ 0'.'305 pixel scale and a 10'5 x 
10'. 5 field of view. Observation were performed over the 
course of 9 observing runs from January 2003 to April 
2006. 

The mean airmass of observation varied between 1.21 
and 1.43 for different field/filter combinations, and ex- 
posures were rarely taken at airmass > 1.6. Standard 
stars from Persson et al. (1998) were observed 2-4 times 
per night except when the conditions were poor. The 
range of airmass values for standard star observations 
was similar to that of science observations. 

The background emission in the NIR is bright, 
non-uniform across the field, and can vary on short 
timescales. Accurate background subtraction requires 
that the telescope be dithered between exposures (see 
§4). Because the brightest objects in an exposure can 
leave residual images in subsequent exposures, using a 
non-regular dither pattern - in which the telescope does 
not repeatedly trace the same sequence of dither posi- 
tions - facilitates removal of artifacts during the reduc- 
tion process. We used an algorithm to generate semi- 
random dither patterns in which the distance between 
subsequent dither positions is maximized. The size of 
the dither box is 45", which is sufficiently large to ob- 
tain good background subtraction in the regions around 
all but the brightest / most extended sources in our fields 
without significantly reducing the area with the highest 
exposure times. 

The typical exposure times at each dither position were 
1 x 100s (coadds x individual exposure time) in J, 4 x 
20s in H, and 4 x 15s in K. The total exposure times 
for each field/filter combination, after discarding images 
with poor quality (§4), are given in Table 1. 

The set of filters used with ISPI changed in April 2004. 
The fields HDFS1 and 1030 were completed using the 
original JHK' filters, whereas HDFS2 and 1255 were ob- 
served using the newer JHK S filters. Both sets of filter 
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Fig. 1. — The filter transmission curves, including atmospheric 
transmission, used for the MUSYC NIR imaging. The solid curves 
are for HDFS1 and 1030 and the dashed curves are for HDFS2 and 
1255. 

transmission curves are shown in Figure 1. Conversions 
between the Vega and AB magnitude systems were cal- 
culated using the SED of Vega, and are given, along with 
the effective wavelengths, in Table 2. Note that the shift 
in effective wavelength between the K' and K s filters is 
small, and much less than the filter widths. In the dis- 
cussion that follows we do not distinguish between the 
two filter sets. 

4. DATA REDUCTION AND IMAGE PROPERTIES 
4.1. Data Reduction 

The data were reduced using a combination of stan- 
dard IRAF 14 tasks, modified IRAF tasks, and custom 
tasks. The core of our reduction procedure was the ex- 
ternal IRAF package XDIMSUM. 15 The basic methods 
are similar to those described in detail by Labbe et al. 
(2003), and are outlined below. 

Near-infrared cameras have a non-negligible dark cur- 
rent, so dark images with the appropriate exposure time 
and number of coadds are subtracted from the science 
images. We use domeflats to perform the flat field cor- 
rection. To remove the background emission from the 
domcflat image, and thus to isolate the uniform illumina- 
tion of the dome screen by the lamp, we follow the stan- 
dard procedure of subtracting a "lamp off" image from a 
"lamp on" image. During the observations, special care 
was taken to keep the count level of the "lamp on" im- 
ages at a reasonable level, well below where non-linearity 
becomes significant on the ISPI camera. Domeflats were 
taken nightly or semi-nightly; they were generally sta- 
ble from night to night, but showed significant variations 
between observing runs. 

14 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

15 http://iraf.noao.edu/iraf/ftp/iraf/extern- 
v212/xdimsum020806 



The background emission in each science frame is 
subtracted, and the resulting images combined, in two 
passes. During the first pass, an image of the back- 
ground emission is created for each science image using 
a running median of the dithered sequence of science im- 
ages; e.g. the background for image 10 is the median 
combination of images 6-9 and 11-14. The positions of 
several stars are used to determine the relative shifts be- 
tween background-subtracted images. The images are 
shifted to a common reference using sub-pixel interpo- 
lation and are combined. Objects in the combined im- 
age are detected using a simple thresholding algorithm. 
These masks are shifted back into the frame of individual 
science exposures, and the background-subtraction pro- 
cess is repeated during the second pass; this time objects 
are masked out during the calculation of the running me- 
dian in order to improve the background subtraction. 

We take several steps to improve the quality of the 
final images. We create a mask of bad pixels in each 
image. An initial list of bad pixels is created using 
the flat-field images. We then inspect each background- 
subtracted image individually; images with severe ar- 
tifacts (e.g. disturbed point spread functions) are dis- 
carded, while others with localized artifacts (e.g. satellite 
trails) are masked using a custom procedure. Additional 
bad pixels in each image are identified using a cosmic 
ray detection procedure or are removed with a sigma- 
clipping algorithm during image combining. The ISPI 
array can retain memory of previous exposures in the 
form of persistence images of bright objects. We create 
a second object mask at the end of the first pass reduc- 
tion, in which only the cores of the brightest objects are 
masked; these are used to mask the pixels that contained 
bright objects during the previous exposure. Finally, we 
optimize the signal-to-noise in the seeing disk by using a 
weighted average during the final image combining step. 
The weights are calculated using 
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{scala * rmsi * FWHMi 



(1) 



where scalei is a constant used to scale the signal in 
image i to the common level, and rmsi is the pixel-to- 
pixel rms measured in a blank region of the unsealed 
image. The values scala and FWHMi are determined 
for each image using a set of bright non-saturated stars, 
which are chosen in locations away from dense regions of 
bad pixels. 

The final products of our reduction procedures are a 
combined image, an exposure time map, and an rms map 
which gives an estimate of the noise level at each pixel. 

4.2. Astrometric Correction and Optical Images 

The observations, reduction, and characteristics of the 
first optical MUSYC data are described in Gawiser et al. 
(2006a). Subsequent data were reduced and analyzed us- 
ing similar methods, and will be described by E. Gawiser 
et al. (2007, in preparation). Most of the optical imag- 
ing was obtained using the 8 CCD MOSAIC II camera 
on the Blanco 4m telescope at CTIO. Each image was 
re-sampled to provide a uniform pixel scale and tangent 
plane projection using stars with known positions from 
the USNO-B catalog. The final rms astrometric errors 
arc estimated to be less than 0'.'2 across the entire field. 
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We use the IRAF tasks GEOMAP and GEOTRAN 
to re-sample the NIR images so that they follow the 
same logical and world coordinate system (i.e. x, y 
pixel coordinates and RA, Dec) as a set of trimmed op- 
tical MUSYC images. This process flips the NIR images 
around the x-axis so that North is up and East is left, 
in accordance with standard practice. The pixel scale is 
changed slightly, from Of! 305 to 0"267. We used a 6th or- 
der fit in x and y, allowing for cross-terms, to adequately 
remove the distortions present in the ISPI instrument; 
a similar high-order fit was also found to be necessary 
by ISPI instrument scientists 16 . We verified that the 
pixel resampling does not introduce systematic errors in 
flux for a set of objects distributed across the images. 
The rms error in stellar positions between the NIR and 
optical images is ~ 0'.'06, much smaller than the typical 
aperture size used for photometry. We also compared 
the positions of stars in our K-band images to those in 
the USNO-B catalog directly. Many of the objects are in 
fact extended, and some others are saturated in our im- 
ages; after removing these objects the rms is 0'.'15 — 0'.'2, 
consistent with the uncertainties of individual USNO-B 
stars. 

The final if-band images are shown in Figure 2. The 
stretch has been adjusted to emphasize faint sources and 
the uniformity of the background. Figure 3 shows a color 
composite image of HDFS1, constructed using the RJK 
bands. The effective seeing in the final images is taken 
as the median of the FWHM for a set of ~ 5 — 10 stars, 
and is given in Table 1. 

4.3. Photometric Calibration and Verification 

Observations of standard stars from Pcrsson et al. 
(1998) were performed nightly except when the condi- 
tions were very poor. Each standard star observation 
used a large 5-point dither pattern, and the telescope 
was de-focused to keep the peak count level in the linear 
regime of ISPI. Following Persson et al. (1998), the in- 
strumental magnitude of the star was measured in a 10" 
diameter aperture; we verified that increasing the aper- 
ture size would change the instrumental zeropoints by 
< 1%. The zeropoint rms from a single dither sequence 
of a standard star is typically 0.01 — 0.03, setting an ap- 
proximate scale for the flatness of the images. We then 
calibrated a set of secondary standard stars in the science 
field with the same size aperture, after applying a small 
correction for the different airmass of standard star and 
science observations. It was not possible to determine ac- 
curate airmass coefficients from our observations, so we 
used the coefficients given by Frogel (1998), which were 
nonetheless found to be consistent with our observations. 
Following usual practice in the NIR, we did not account 
for color terms in the photometric calibration. 

There was at least one night for each field/filter com- 
bination where the observing conditions were sufficient 
to determine an accurate calibration of the secondary 
standards. For field/filter combinations where there was 
more than one night of satisfactory calibrations, the 
agreement between nights is typically 0.01 — 0.02. While 
it is difficult to place constraints on possible systematic 
uncertainties, the internal precision of our final zero- 
points is < 0.03 in every case. The zeropoints are given 

16 see http://www.ctio.noao.edu/instruments/ir_instruments/ispi/ 



in Table 3. 

The quality of our photometric calibrations was ver- 
ified in several ways. There is a ~ 10 arcmin 2 overlap 
region between two of our fields, HDFS1 and HDFS2, 
where the exposure time is at least half of the total. As 
these fields were observed on different runs, using differ- 
ent NIR filters, and the photometric zeropoints were de- 
termined independently, comparing the total magnitudes 
(§ 5) of objects in this region from the two fields allows a 
check of the internal consistency of our photometric cali- 
brations. We find that the median offset is less than 0.02 
for bright objects in each of J, H, and K. It should be 
noted that these objects lie near (opposite) edges of the 
ISPI detector in the two fields, where the flat-fielding er- 
rors may be larger; the good agreement between the two 
fields is encouraging. We do not attempt to adjust the 
zeropoints to obtain better agreement. 

The MUSYC HDFS1 field contains the - 4.5 arcmin 2 
FIRES HDF-S field (Labbe et al. 2003), allowing for an 
additional check of our photometry The FWHM of the 
FIRES image is ~ 0'.'45. Rather than using the total 
magnitudes from the FIRES catalog, we convolved the 
public FIRES image to match the MUSYC point spread 
function (PSF), and directly compared the aperture pho- 
tometry for a set of compact objects. We found the agree- 
ment to be better than 0.02 magnitudes in each of JHK. 

A final check comes from comparing our total mag- 
nitudes to the aperture-corrected magnitudes from the 
public 2MASS point-source catalog (Skrutskie et al. 
2006). We selected a set of matching objects from the 
three MUSYC deep NIR fields, rejecting objects with 
low S/N in the 2MASS photometry, that were affected 
by non-linearity/saturation in the MUSYC photometry, 
or were blended with other objects. The mean offset 
ranges between 0.01-0.04 for HDFS1 and HDFS2, gener- 
ally consistent with the errors. However, adopting zero- 
points based on the 2MASS photometry would decrease 
the level of agreement for objects in the overlapping re- 
gion of HDFS1 and HDFS2 for two out of three NIR fil- 
ters, and would decrease the agreement between HDFS1 
and the FIRES HDF-S photometry for all three filters. 
For these two fields, we conclude that our flux calibration 
and the quoted zeropoint uncertainties are reliable, and 
are better than could be obtained by calibrating directly 
from the 2MASS catalogs. 

However, the agreement between 2MASS photome- 
try and MUSYC photometry is worse for the 1030 and 
1255 fields. The mean difference between 2MASS and 
MUSYC 1030 photometry is a remarkably consistent 
~ 0.050 ±0.015 in each of JHK, where the uncertainties 
are given as the standard deviation of the mean. The 
sense of the disagreement is that the stars are brighter 
in the 2MASS catalog than in the MUSYC catalog. The 
disagreement is of similar size, but in the opposite di- 
rection, for 1255. This is worrisome, as the level of dis- 
agreement for these fields in all three filters is worse than 
in any of the filters for HDFS1 and HDFS2, and more- 
over the differences are systematic across the filters. We 
have verified that this disagreement is not caused by ob- 
vious mistakes in our photometric calibrations. The sys- 
tematic nature of the offsets may indicate that they are 
caused by aperture corrections in 2MASS, rather than ze- 
ropoint errors. Because our calibrations appear to work 
very well for HDFS1 and HDFS2, we do not adjust our 
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Fig. 2. — The deep K-band images from the MUSYC. At top, the 
a single 19'5 X 10'3 image. The other fields are ~ 10.'3 on a side. 

zeropoints using 2MASS photometry; but without other 
data to compare with, we cannot be certain that our 
< 0.03 magnitude zeropoint uncertainties for these fields 
are reliable. 

4.4. Noise Properties and Limiting Depths 

The flux uncertainty within an aperture has a con- 
tribution from the photon statistics from astronomical 
objects, as well as a contribution from background noise, 
which is due to sky emission, read noise, etc. The stan- 
dard method of determining the contribution of back- 
ground noise, by e.g. SExtractor (Bcrtin & Arnouts 
1996) or the APPHOT package of IRAF, is to scale 
the rms flux value of background pixels by the square 
root of the number of pixels within the aperture, N p i X . 
This scaling is appropriate if adjacent background pix- 



individual HDFS1 and HDFS2 images have been combined to create 

els are uncorrelated, however correlations will be intro- 
duced by imperfect subtraction of background emission, 
extended wings from bright objects, undetected sources, 
re-pixelization during the reduction procedure, and ar- 
tifacts in the images. In the limiting case of perfect 
correlations between the pixels within an aperture, the 
background noise aback will scale as N pix . Thus it might 
be expected that the true scaling will be a oa ck oc N^ ix , 
where 0.5 < (3 < 1. 

We follow Labbe et al. (2003) in characterizing the 
noise properties of our images by summing the counts 
in apertures distributed randomly over empty regions 
of each image. The rationale and method are discussed 
more fully in Gawiser ct al. (2006a). Figure 4 shows the 
distribution of fluxes in l'/5 and 2" diameter apertures. 
These distributions are well-approximated by Gaussians. 
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Fig. 3.— An RJK composite of HDFS1 



Figure 5 shows how the width of the best-fitting Gaussian 
changes with aperture size, along with predictions from 
the -J ' Npi x and N p i x scalings that are described above. 
It is apparent that neglecting correlations between pixels 
would cause a significant underestimation of background 
fluctuations. A single power law provides a good fit to 
the relationship between N P i X and aback for the aperture 
sizes of interest; the typical power law index is (5 ~ 0.6. 

We estimate the flux uncertainties for each object in 
each band as 

F 

° 2 = a bach + ^jf^, (2) 

where F is the flux in ADU, GAIN is the total effec- 
tive gain, and aback is estimated for an aperture with the 
appropriate size. In the case of the elliptical Kron aper- 
tures (see § 5) we estimate aback for a circular aperture 
with the same area. 



The characteristics of the final images can be found 
in Table 3. This table includes an estimate of the 5a 
point source limiting depths. These values are calculated 
using the background fluctuations in the color aperture 
(see §5.3), and with a 0.75 magnitude aperture correction 
applied to account for the ~ 50% point-source flux that 
falls outside this aperture. An alternate analysis of the 
limiting depths can be found in § 5.5. 

5. SOURCE DETECTION AND PHOTOMETRY 

5.1. Source Detection 

We used the SExtractor v2.4.3 software (Bertin & 
Arnouts 1996) to detect objects in the if -band im- 
ages. The two parameters that affect the sensitivity 
of source detection are DETECT.THRESH, which was 
used to specify the detection threshold in units of the 
mis of background pixels, and DETECT_MINAREA, 
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Fig. 4. — Histograms of background fluctuations in 1'.'5 and 
2" diameter apertures in the HDFS1 A"-band image. The larger 
aperture has the broader distribution of fluxes. The dashed curves 
are the best-fitting Gaussians. 
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Fig. 5. — The rms value of background fluctuations within an 
aperture as a function of the aperture size for HDFS1 K. The 
filled circles illustrate the measured values, while the solid curve 
is a power law fit. The dashed curves show the expected scaling 
from the measured pixel-to-pixel rms relation in the case of no pixel 
correlations (lower) and perfect correlation of all pixels within each 
aperture (upper). 

which specifies the number of adjacent pixels that must 
meet this threshold. Additionally, SExtractor option- 
ally filters the detection image with a convolution ker- 
nel prior to detection in order to enhance the detec- 
tion of faint objects. Rather than use detailed simula- 
tions to find the optimum set of parameters for SEx- 
tractor, wc made use of the ultradeep K-band image 
of HDFS from the public FIRES survey (Labbe et al. 
2003). The difference in depth between the MUSYC 
HDFS1 if-band image and the FIRES HDF-S image 



(K ~ 21.0 and K ~ 24.3, respectively) allows for a 
clear determination of which objects detected by SEx- 
tractor in HDFS1 are real and which are noise peaks. 
We increased the DETECT.THRESH parameter slightly 
from the default value in order to eliminate all spurious 
sources in the FIRES region, but without significantly 
affecting the completeness of our catalogs. The final de- 
tection parameters were DETECT_MINAREA=5 pixels 
and DETECT.THRESH = 1.6, and we used a Gaussian 
filter corresponding to the PSF of the detection image. 
We detected objects in rms-normalized if-band images, 
which were created by dividing the i-T-band science im- 
ages by the rms maps (see §4). Detection in normal- 
ized images is important because the dithering procedure 
used during observations means that the outskirts of the 
science images are noisier, which would otherwise result 
in a large number of spurious sources. 

5.2. PSF matching 

In order to measure accurate colors, all optical and NIR 
images were PSF matched using Gaussian convolution 
kernels. This is necessary to ensure that the aperture 
used for color measurements (§5.3) contains a constant 
fraction of the total light for every object in every band. 
The systematic effects on colors that are introduced by 
imperfect PSF matching can be estimated by dividing 
the stellar growth curves, which were created for each 
image using a set of bright non-saturated stars evenly 
distributed over the images. We optimized the width 
of the convolution kernels by minimizing the difference 
in growth curves; we found that simple Gaussians are 
capable of reducing the systematic errors to the ~ 1 — 
2% level for apertures as large or larger than the color 
apertures. The exception is the z-band, which has a PSF 
shape that is sufficiently different from the other bands 
that a Gaussian kernel yields errors that are ~ 2 — 3%. 

Most of the images have FWHM from ~ 0'/8 - l'.'O. 
However the U and B bands are slightly worse, with 
values as high as l'/4. Degrading all the images of a 
given field to the broadest PSF would reduce the S/N 
significantly in every band. For the images with narrow 
PSFs, we smooth to ~ 1" in order to measure uniform 
colors. Images with broader PSFs are treated differently. 
We measure e.g. the U — K color using a version of the 
K image that is degraded to the same PSF as U, and 
add this color to the AT-band flux in the color aperture 
of the ~ 1" K image. It is important to note that this 
procedure may not be appropriate for extended objects 
with strong color gradients. This is because the U — 
K color is measured for a larger region of the galaxy 
than is e.g. the J — K color, with the effect that the 
inferred U — J color will be different than what could be 
measured in any aperture. This effect should not be a 
significant concern here, because distant galaxies are at 
best marginally resolved in our ~ 1" seeing and observed 
color gradients are minimal. Indeed, we have verified 
that any errors introduced by this procedure are at no 
more than the percent level by directly measuring the 
U — J color (using a smoothed version of the J image) 
and comparing it to what is inferred from the U — K and 
J — K colors. 

5.3. Photometry and Colors 
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We use SExtractor in dual-image mode to detect ob- 
jects in K and to perform photometry in all bands. In 
this subsection we describe our methods to determine 
"total" flux in K, and to measure high signal-to-noise 
(S/N) colors. The total flux in any other band can be 
calculated directly from these quantities. 

We estimate the total if -band flux of an object from 
the SExtractor AUTO photometry, which uses a flexible 
elliptical Kron-like aperture. Although this aperture con- 
tains most of the flux for bright objects, in practice the 
aperture can become quite small for faint sources-even 
for the same light profile-so a larger fraction of light may 
be missed. Therefore we convert the AUTO flux to total 
flux by applying an aperture correction. The aperture 
correction is calculated using the median stellar growth 
curve of a set of bright stars, for a circular aperture with 
the same area as the AUTO aperture. We note that 
this procedure may introduce a mild bias in the total 
flux of faint, extended sources, because their light pro- 
files may not follow that of a point source outside the 
AUTO aperture; however these sources will have highly 
uncertain flux measurements no matter what procedure 
is used. The aperture correction for the AUTO aper- 
ture can reach ~ 0.5 magnitudes for some of the faintest 
sources. 

Because accurate color measurements are necessary for 
photometric redshift calculations and modeling of stel- 
lar populations, we wish to optimize the S/N. Since the 
noise is a strong function of aperture size (Fig. 5), we 
measure the colors in smaller apertures than the AUTO 
aperture that is used to estimate the total flux. For a 
point source with a Gaussian PSF, and for uncorrelated 
background noise, the aperture with optimal S/N has 
diameter 1.35 x FWHM (Gawiser et al. 2006a). In real- 
istic situations there are competing effects which change 
the size of the optimal aperture. Relative to the ideal- 
ized case, the broader wings of our PSFs would suggest 
a larger aperture to obtain the same signal, while the 
noise correlations suggest a smaller aperture to keep the 
same level of noise. The optimal S/N aperture for a point 
source can be found by dividing the stellar growth curve 
by the aback curve shown in Figure 5. We find that the 
optimal aperture typically has diameter ~ 1.1 — 1.4 times 
the FWHM, depending on the filter (see also Gawiser et 
al. 2006a). Using a very small aperture for color deter- 
minations presents several problems. Accurate colors re- 
quire that a similar fraction of the flux from an object is 
contained within the color aperture in each filter, which 
is easier to achieve with larger apertures. This is because 
a larger aperture contains a larger fraction of the flux, 
leading to a smaller relative difference in the fraction of 
flux. Secondly, one of the primary scientific goals of the 
deep NIR MUSYC survey is the study of high redshift 
galaxies, which may not appear as pure point-sources in 
our ~ 1" images. Larger apertures also reduce the effects 
of the small variations in PSF across the images, and the 
residual geometric distortions at the edges of the images. 
We choose, as a compromise, apertures with diameter 
~ 1.5 times the stellar FWHM. These apertures contain 
~ 50% of the light from a point source. We verified that 
the S/N in these apertures is > 95% of the S/N in the 
optimal aperture for every image. 

5.4. Catalog format 



The photometry in the if -selected catalogs is presented 
in units of flux, normalized so that the zeropoint is 25 on 
the AB system. The use of flux, rather than magnitudes, 
avoids the problem of converting the measured flux un- 
certainties into magnitude uncertainties, the problem of 
asymmetric magnitude uncertainties for low S/N objects, 
and the loss of information for objects that have nega- 
tive measured fluxes. The photometry is not corrected 
for Galactic extinction. 

Object detection and measurement of geometrical 
parameters are performed by SExtractor in the rms- 
normalized if image. The adjacent HDFS1 and HDFS2 
fields were treated separately, but a list of objects that 
appears in both catalogs is available. Version 3.1 of 
the if -selected catalogs are given in the following format: 

Column 1: SExtractor ID number, starting with 1 

Columns 2-3: x and y barycenters 

Columns 4~5: RA and DEC (decimal degrees; J2000) 

Column 6: internal MUSYC field code (4=HDFS1, 
5=HDFS2, 6=1030, 7=1255) 

Columns 7-24: Flux density and uncertainty in the 
color aperture, in the order U BVRIzJHK 

Columns 25-26: if-band flux density and uncertainty 
in the "total" aperture 

Columns 27-35: Exposure time weight in bands 
U BVRIzJHK, normalized to the weight of the median 
object 

Columns 36-37: Diameters of color and AUTO aper- 
tures. The diameter of the AUTO aperture is taken as 
the geometric mean of the major and minor axes (arcsec) 

Columns 38: SExtractor blending flag 2 - object 
was originally blended with another object (l=blended, 
0=unblended) 

Columns 39: SExtractor blending flag 1 - object was 
blended with another object strongly enough to sig- 
nificantly bias AUTO photometry (l=blended, 0=un- 
blcnded) 

Column 40: Half-light radius (arcsec) 
Column 41: Ellipticity 

Column 4% : Position angle (degrees), measured 
counter-clockwise from North 

Column 43: Aperture correction to convert AUTO flux 
to total flux 

Column 44 : SExtractor CLASS_STAR parameter 

Column 45: Maximum SExtractor flag, for bands 
UBVRIzJH 

Column 46: SExtractor flag, if -band 

The catalogs are available in the electronic edition of 
the journal article. 

5.5. Completeness 

We estimate the completeness of our catalogs as a func- 
tion of magnitude by attempting to detect simulated 
point sources. The point sources are created by extract- 
ing a bright, non-saturated star, scaling it to the desired 
flux level, and inserting it at random locations in the 
central well-exposed regions of our if-band images. We 
then attempt to detect the stars using the same SExtrac- 
tor settings described in § 5. Figure 6 shows the resulting 
completeness curves as a function of magnitude. 

Because the simulated point sources are placed at ran- 
dom locations, some fraction will fall on or near enough 
to other objects that SExtractor does not properly de- 
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Fig. 6. — The X-band completeness curves. Point sources 
were inserted at random locations in the central region of the four 
MUSYC fields, where the exposure time is > 95% of the total. The 
completeness is defined as the fraction of simulated sources that 
were recovered, as a function of total magnitude of the simulated 
source. The completeness values arc higher if the simulated point 
sources are inserted in empty regions of the images (see Table 4). 

blend them. This effect causes the plateau in the com- 
pleteness curves at brighter magnitudes to have some 
slope; evidently, this problem is worse at fainter mag- 
nitudes. We repeated the simulations, this time insert- 
ing point sources in locations that avoid other (real or 
simulated) objects. The completeness is now equal to 
unity to at least K = 20.5 in all fields. Table 4 shows 
the 90% and 50% completeness limits for the deep K- 
selected MUSYC catalogs for both the "unmasked" and 
"masked" simulations. 

We note that the completeness limits from the un- 
masked and masked simulations are similar, reflecting 
the uncrowded nature of the fields to our depth in K. 
The masked simulations may provide a better estimate 
of the sensitivity of our images, while the unmasked 
(shallower) simulations should be used when assessing 
the completeness for actual astronomical objects because 
such objects do not avoid each other on the sky. We also 
note that the completeness for extended objects will be 
lower than for the simulated point sources, and would in 
principle be a function of inclination, morphology, and 
size. 

6. NUMBER COUNTS 

Figure 7 shows the surface density of objects as a func- 
tion of magnitude, excluding objects classified as stars 
(§7), for each of our four fields. No completeness correc- 
tions have been applied. The number counts are calcu- 
lated in 0.5 magnitude bins using only the image area 
with > 95% of the total X-band exposure time. The 
errorbars assume Poisson statistics, which should un- 
derestimate the true uncertainties because objects arc 
clustered. The four fields are generally consistent over 
18 < K < 20, where the S/N is high and differences in 
depth do not affect the number counts. It is noteworthy 
that the 1030 field has the highest density of objects in 



this magnitude range; this field also has the highest den- 
sity of galaxies with M > lO 11 ^© at 2 < z p hot < 3 of any 
of the fields studied by van Dokkum et al. (2006), sug- 
gesting the possibility of significant galaxy over-densities 
in this redshift range. 

The upper panel of Figure 8 shows the average num- 
ber counts of the four fields, along with counts drawn 
from the literature. The completeness corrections de- 
scribed above have been applied to the MUSYC points. 
Note that these completeness values provide a simplis- 
tic correction when dealing with number counts; a more 
sophisticated correction would account for the difference 
between the measured and intrinsic magnitudes of the 
artificial sources. This can be a significant effect be- 
cause the rising slope of number counts means that more 
sources would scatter to brighter magnitudes than to 
fainter magnitudes due to noise fluctuations. Further- 
more, there may be biases in the photometry of the 
faintest sources (see also Forster Schreiber et al. 2006). 
Neither do we attempt to correct for spurious detections, 
which would only be significant at the faintest magni- 
tudes, or estimate completeness correction for extended 
sources. For these reasons we only extend the average 
number counts to the bin centered at K = 21, where the 
completeness correction begins to become significant. 

The best fit logarithmic slope d{logN)/dM to the 
MUSYC number counts isa« 0.31 over 18 < K < 20. 
To illustrate deviations from this power law in different 
magnitude ranges, in the lower panel of Figure 8 we have 
divided the observed number counts by this fit. The de- 
viation from unity in this panel at brighter magnitudes 
illustrates a change in the d(logN)/dM relation. How- 
ever, we do not find evidence for a sharp break in the 
power law slope at K ~ 17.5, as has been reported by 
Cristobal- Ilornillos et al. (2003). These authors interpret 
the galaxy number counts in terms of models of galaxy 
evolution, and suggest that the break can be reproduced 
by models with late star formation in massive galaxies ( 
z < 2; see also Eliche-Moral et al. 2006). There is also 
some evidence from deeper surveys that the slope flattens 
at K > 21; for instance, the data of Forster Schreiber et 
al. (2006) indicate a w 0.20 over 21 < K < 23. 

7. PHOTOMETRIC REDSHIFTS AND STAR 
CLASSIFICATION 

We calculate photometric redshifts using the methods 
described by Rudnick et al. (2001, 2003). Briefly, non- 
negative linear combinations of galaxy templates are fit 
to the observed spectral energy distributions. The tem- 
plates include the four empirical templates of Coleman, 
Wu, & Weedman (1980) and the two empirical star- 
burst templates of Kinney et al. (1996), all of which have 
been extended into the UV and NIR using models. We 
also include lOMyr and lGyr old single stellar population 
templates generated with the Bruzual & Chariot (2003) 
models, as the empirical templates are derived from low- 
rcdshift galaxies and do not adequately describe some of 
the high- redshift galaxies in the MUSYC catalogs. We 
do not allow for additional reddening in the models. The 
photometric redshift uncertainties are calculated using 
Monte Carlo simulations, in which the observed fluxes 
are varied within the photometric uncertainties. We cor- 
rected for Galactic extinction when calculating photo- 
metric redshifts (Schlegel et al. 1998). 
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Fig. 7. — The raw number counts in the four MUSYC fields, ex- 
cluding objects classified as stars. Poisson errorbars are only shown 
for the HDFS1 field; the other fields have comparable uncertainties. 
No completeness correction has been made. 
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Fig. 8. — Top: The average galaxy number counts from the four 
MUSYC fields, along with a compilation of results from the litera- 
ture. The MUSYC points have been corrected for incompleteness, 
but only points where this correction is small are shown. Bottom: 
The number counts have been divided by a power law with index 
a Ri 0.31, which provides the best fit to the MUSYC points over 
K < 18 < 20. The procedure highlights the increasing logarith- 
mic slope of the number counts at K < 18. The slope may also 
decrease slightly at K > 21. 



Fig. 9. — Comparison of photometric and spectroscopic rcdshifts. 
Errorbars represent the 68% confidence intervals deteremined from 
Monte Carlo simulations. The median |Az|/(l + z) is 0.044. As 
described in the text the spectroscopic redshifts for six galaxies 
are shown with some uncertainty, as they are derived from fits 
to the stellar continuum rather than direct observations of emis- 
sion/absorption lines. 

We have compiled a list of spectroscopic redshifts from 
the ongoing MUSYC spectroscopic program (P. Lira et 
al. 2007, in preparation; Kriek et al. 2006b, c) and the 
NASA Extragalactic Database. Figure 9 compares the 
spectroscopic and photometric redshifts. The six spec- 
troscopic redshifts taken from Kriek et al. (2006b) are 
shown with some uncertainty, as they are derived from 
fits to the stellar continuum rather than direct measure- 
ments of emission lines. These galaxies have prominent 
Balmer/4000A breaks that strongly constrain the fits, so 
there is little chance of catastrophic failure in the redshift 
estimates. In all there are 130 spectroscopic redshifts, 35 
of which are at z > 1.5. 

> 



The median of |Az|/(l + z) is 0.044 (0.07 at 



'spec 



1.5). Catastrophic outliers, which we define as those ob- 
jects with | Az|/(l+z) greater than five times the median, 
comprise ~ 8% of the sample. 

We find that the SExtractor CLASS_STAR parame- 
ter does not give a reliable stellar classification for all 
objects for our i\-selected catalogs; some objects with 
a high stcllarity parameter show obvious extended pro- 
files in the optical bands. Rather than relying on SEx- 
tractor to identify stars, we fit the NextGen stellar at- 
mosphere models (Hauschildt et al. 1999) to all objects 
in the MUSYC catalogs. Stars are identified as those 
objects that are better fit by stellar atmosphere models 
than by galaxy templates, and which do not have clearly 
extended profiles in either the NIR or optical bands. 

8. ANALYSIS OF COLOR SELECTION TECHNIQUES 

In this section we use the MUSYC data and photomet- 
ric redshifts to shed light on the relation between galaxies 
that are selected with the BM/BX/LBG criteria (Steidel 
et al. 2003; Adelberger et al. 2004), the DRG criterion 
(Franx et al. 2003), and the BzK criteria (Daddi et al. 
2004) . These selection techniques were designed with dif- 
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ferent types of galaxies in mind, and for different sets of 
observational constraints, so it is not immediately clear 
how they should relate to each other. First, we give a 
brief discussion of these techniques and describe how we 
implement them. Then we show the locations of galaxies 
selected using these criteria on the BzK and J — K diag- 
nostic diagrams, as has been done previously by Reddy 
et al. (2005) for a smaller spectroscopic sample. Finally, 
we present photometric redshift distributions and rest- 
frame optical colors for these these galaxies. 

The galaxies discussed in this section are drawn from 
the four deep MUSYC fields. We select galaxies with 
K < 21 and require a minimum if-band exposure time 
weight of 0.6. 

8.1. The color selection criteria 

Lyman break and BM/BX galaxies: The classical "£/- 
dropout" technique has proven very effective at iden- 
tifying the so-called Lyman break galaxies (LBGs) at 
z ~ 3. Adelberger et al. (2004) introduced the BM 
and BX selection criteria, which arc designed to select 
galaxies with similar SEDs to the LBGs, but at red- 
shifts z ~ 1.7 and z ~ 2.3, respectively. These galax- 
ies are selected to be blue in the rest-frame UV. Stei- 
del and collaborators also apply an TZab < 25.5 limit 
to ensure a reasonable degree of photometric and spec- 
troscopic completeness. Because these galaxies are com- 
paratively bright in the optical and frequently have de- 
tectable emission/absorption lines, they are well-suited 
to spectroscopic follow-up (e.g. Steidel et al. 2003). How- 
ever, the BM/BX/LBG techniques may miss very dusty 
star-forming galaxies or galaxies with little ongoing star 
formation. 

The specific color criteria used by Stcidcl et al. (2003) 
and Adelberger et al. (2004) are based on U n GlZ col- 
ors, but these filters were not used as part of MUSYC. 
Rather than develop analogous criteria using our filter 
set, we calculated synthetic U n GTZ colors from the best- 
fitting models (see also van Dokkum et al. 2004; Daddi 
et al. 2004; van Dokkum et al. 2006). These models were 
generated using the Bruzual & Chariot (2003) popula- 
tion synthesis code. In order to provide the best fits we 
allowed for a range of star formation histories and ex- 
tinction values, but the redshifts were fixed at the values 
determined in §7. 

In the discussion that follows, we draw the distinction 
between BM/BX/LBGs that meet the ft A B < 25.5 limit 
and those that do not. We also recall that these galax- 
ies are usually selected in the optical rather than in K, 
so our conclusions may not directly apply to the wider 
population of BM/BX/LBGs . 

Distant red galaxies: The observed J — K color probes 
the rest-frame optical at z ~ 2 — 4. Stars and galaxies at 
lower redshifts tend to have blue J — K colors, so the sin- 
gle J — K > 2.3 criterion for distant red galaxies (DRGs) 
should primarily select red galaxies at z > 2 (Franx et 
al. 2003; van Dokkum et al. 2003). DRGs are selected in 
K , rather than in the optical, which is closer to a mass 
selection than a selection on unobscured star formation. 
Because DRGs are not expected to include galaxies with 
blue colors or low stellar masses, they may be comple- 
mentary to BM/BX/LBGs. The relative depths of the 
MUSYC JK images are well-matched for the identifica- 
tion of these red galaxies. 



BzK galaxies: Daddi et al. (2004) use two sets of color 
criteria involving the B, z, and K bands to select star 
forming galaxies at 1.4 < z < 2.5 (hereafter, sBzK galax- 
ies) and passive galaxies over the same redshift range 
(hereafter, pBzK galaxies). Together, these selection 
techniques are designed to identify a nearly complete 
sample of galaxies in the targeted redshift range. The 
BzK criteria were designed to be used on if-selected ob- 
jects, so they may not include low-mass galaxies, even 
when bright in the rest-frame UV. 

The B, z, and K s filters used by Daddi et al. (2004) to 
define the BzK color cuts are very similar to the filters 
used by MUSYC, with the most significant difference be- 
ing that the MUSYC K filters extend slightly further to 
the blue; this introduces a small 300A shift in effective 
wavelength. We make no effort to account for differences 
between the filter sets. Many of the redder galaxies are 
very faint or undetected in either B or z; for galaxies 
that are < 3<r detections, we use synthetic magnitudes 
that are calculated directly from the best-fitting galaxy 
templates. The most significant effect of this procedure 
is to place objects that are undetected in either B or z 
in a reasonable location on the BzK diagram. 

The use of synthetic B, z magnitudes for some galaxies 
is necessitated by our limited depth. In fact, the syn- 
thetic B magnitudes suggest that many of our reddest 
-ft'-selected galaxies are so faint in B that they could not 
be observed in even the deepest astronomical images. 
The same is not true for z: the synthetic z magnitudes 
suggest that the reddest galaxies in our catalogs would 
still be detected at > 2a in the study recently performed 
by Kong et al. (2006), which was optimized for the use 
of the BzK selection technique. We note that, to reach 
this depth, those authors observed with an 8m telescope. 
It follows that, if our B and z images reached the depth 
acheived by Kong et al. (2006), a significant number of 
galaxies would have lower limits in B — z. As this color 
is necessary to distinguish between sBzKs and pBzKs 
(see Fig. 10), it may not always be possible to make this 
distinction even with very deep optical images; it will, 
however, be clear that the galaxies are either sBzKs or 
pBzKs. For the sake of simplicity in the discussion that 
follows, we assume that the distinction between sBzKs 
and pBzKs can always be made. 

8.2. Observed Colors 

8.2.1. The BzK color-color diagram 

Figure 10 shows the BzK color-color diagram, marking 
the regions that separate the pBzKs, sBzKs, and stars. 
The left panel highlights the location of BM/BX/LBGs 
on this diagram. The majority of these galaxies fall 
within the sBzK region, with ~ 20% falling outside 
(Reddy et al. 2005 find a similar result). This is par- 
tially due to photometric errors, which can be of signif- 
icant concern for galaxies lying near the bottom of the 
sBzK selection region. Since the sBzK selection window 
was designed to isolate galaxies at z > 1.4, it is not a sur- 
prise that our photometric redshifts suggest that ~ 90% 
of the BM/BX/LBGs that fall outside this window are 
at z < 1.4. Roughly half lie at 1 < z < 1.4 and half 
are the BM/BX/LBG interlopers at z < 1. Additionally, 
some of the BM/BX/LBGs lie in the region of the BzK 
color-color diagram that is supposed to isolate stars (see 
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Fig. 10. — The BzK color-color diagram used by Daddi et al. (2004) to separate star forming galaxies at z > 1.4, passive galaxies at 
2 > 1.4, and stars. Left: The filled blue circles are selected using the BM, BX, or LBG criteria, while the open purple circles are galaxies 
that meet these color criteria but are fainter than the typical limit of TZ\b < 25.5. Small black points are the remaining MUSYC galaxies. 
Right: The filled red circles are DRGs. The larger filled green circles are the nearly passive galaxies from Kriek et al. (2006b). The black 
stars mark the objects identified as stars using our SED-fitting and morphological criteria. 



also Reddy et al. 2005, Fig. 12). 

While the sBzK selection appears to identify most of 
the BM/BX/LBGs at z > 1.4, these optical selection 
techniques are less efficient at identifying sBzKs: roughly 
35% of the sBzKs at our limit of K < 21 are not selected 
by any of the optical criteria. This quantity increases to 
~ 60% at K < 20. These values are in rough agreement 
with the results from the smaller spectroscopic samples 
of Daddi et al. (2004) (see their Fig. 13) and Reddy et 
al. (2005), and reflect that fact that a subset of sBzKs 
have red rest-frame UV colors. Finally, there is virtually 
no overlap between the BM/BX/LBGs and pBzKs. 

The right panel of Figure 10 shows the DRGs on the 
BzK diagram. Most of the DRGs lie in either the sBzK 
region or pBzK region. We note that many of the DRGs 
are very faint in the optical, so very deep B and z imaging 
would be necessary to accurately select these galaxies 
using the BzK techniques (see 8.1). Even so, as shown 
by Franx et al. (2003) many red galaxies escape detection 
in even the deepest optical images. 

The right panel of Figure 10 also shows the locations 
of 9 galaxies with strongly suppressed star formation, 
taken from the sample of Kriek et al. (2006b). The evi- 
dence for low ongoing star formation comes from the lack 
of detectable rest-frame optical emission lines in deep 
NIR spectroscopy, and from stellar population synthesis 
modeling which indicates star formation rates of order 
~ lM^/yr. Despite their nearly passive nature, Fig- 
ure 10 shows that only 3 of these 9 galaxies are classified 
as pBzKs-with another galaxy very close to the pBzK se- 
lection window-and the remaining galaxies are sBzKs. It 
is possible that this reflects a redshift-dependence in the 
effectiveness of the pBzK criterion, even within the range 
of redshifts for which this criterion was designed: of the 
9 nearly passive galaxies that lie within or very near to 
the pBzK region, 4 are at 2.4 < z < 2.6, whereas the 
others are at 2.0 < z < 2.4. We also note that the stellar 
population modeling performed by Kriek et al. (2006b) 
suggests that the former 4 galaxies may have the oldest 



relative ages. They find age/r > 25 for these 4 galaxies, 
where r is the star formation e-folding timescale, and the 
remaining 5 have age/r < 25. Thus it could also be that 
pBzK criterion only works well for the oldest subset of 
nearly passive galaxies. 

8.2.2. The J — K versus K color-magnitude diagram 

Figure 11 shows the J — K versus K diagram, along 
with a line at J — K = 2.3 which marks the limit 
for DRGs. The left panel shows that only ~ 10% 
of BM/BX/LBGs are classified as DRGs. Conversely, 
~ 32% of DRGs satisfy one of the BM/BX/LBG criteria. 
Of the galaxies that meet both sets of criteria, ~ 55% are 
fainter than the IZab = 25.5 limit used in many ground- 
based samples of optically-selected galaxies. 17 Thus the 
overlap between DRGs and BM/BX/LBGs is small, and 
has a significant contribution from galaxies that are too 
faint to be included in typical optical surveys (i.e. they 
have 7^ab > 25.5) (see also Reddy et al. 2005; van 
Dokkum et al. 2006). 

The right panel shows that sBzKs and pBzKs span a 
range in J — K color. This panel also shows the loca- 
tion of the nearly passive galaxies described by Kriek et 
al. (2006b). The DRG criterion identifies 6/9 of these 
galaxies. This incompleteness may be due to photomet- 
ric uncertainties, as 2/3 of the remaining galaxies are 
within ~ lcr of the DRG limit. As shown above, only 
3/9 are selected by the pBzK criterion. 

8.3. Photometric redshift distribution 

Figure 12 shows the photometric redshift distributions 
for each of the color criteria discussed above. The dis- 
tribution of DRGs shows a prominent gap in the region 
z ~ 1.4 — 1.8. It is not clear whether this is an artifact 
in our photometric redshift or a real feature. Defining 
interlopers as those objects at z < 1.8, the interloper 



17 These results do not change significantly if the z 
galaxies are removed from the analysis. 
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Fig. 11. — J — K versus K. The lines at J — K = 2.3 illustrate the DRG limit. Left: The filled blue circles are selected using the 
BM, BX, or LBG criteria, while the open purple circles are galaxies that meet these color criteria but are fainter than the typical limit 
of "R-ab = 25.5. Small black points are the remaining galaxies. A slight majority of BM/BX/LBGs which are also classified as DRGs are 
fainter than this limit. Right: Filled blue circles are sBzKs, and open red circles are pBzKs. The larger filled green circles are the nearly 
passive galaxies from Kriek et al. (2006b). Two of the three passive galaxies that are not classified as DRGs are within ~ la of the DRG 
limit. The black stars mark the objects identified as stars using our SED-fitting and morphological criteria. 




z phot 

Fig. 12. — The photometric redshift distribution of K < 21 
galaxies selected according to different selection criteria. The filled 
histograms are for galaxies with K < 20. The dashed, dotted, 
and dot-dashed histograms on the bottom panel are for the BM, 
BX, and LBG criteria, respectively. The upper dashed lines in 
the lower panels show the redshift distribution of all objects in the 
MUSYC catalogs, illustrating that the sBzK and BM/BX/LBG 
criteria select most of the galaxies in the higher redshift bins. 

fraction is ~ 15% at K < 21, ~ 20% at K < 20, and 
~ 50% at K < 19. 

The photometric redshift distribution of pBzKs is 
shown in the second panel of Figure 12. The photo- 



metric redshifts indicate that this selection technique is 
indeed effective at isolating galaxies at z > 1.4. The in- 
terlopers lie primarily at 1.2 < z < 1.4, with very few 
at lower redshift. As shown in the third panel, the sBzK 
galaxies have an approximately similar redshift distri- 
bution but with a substantial number reaching down to 
z ~ 1 — 1.2. The relatively minor differences in photomet- 
ric redshift distributions suggest that these techniques 
may be used to select complementary samples of galax- 
ies, although the difference in the number of interlopers 
should be taken into account: 10% of the pBzKs and 
25% of the sBzKs lie at z < 1.4 (these numbers increase 
to 14% and 30%, respectively, for K < 20, and to 22% 
and 30% for K < 19). Also, as shown above, galaxies 
with quenched star formation may be classified as sBzKs. 
For these reasons, it appears that precise comparisons of 
the global properties of star forming and passive galaxies 
(e.g. number density, stellar mass, or clustering) using 
only BzK photometry may be of limited usefulness (see 
Kong et al. 2006). 

While the BzK criteria were designed to select galax- 
ies at 1.4 < z < 2.5, our photometric redshifts suggest 
that they identify a fairly complete sample of galaxies 
at larger redshifts, up to z ~ 3.5. Together, these crite- 
ria are very effective at selecting high redshift galaxies: 
~ 93% of the galaxies in our catalogs at 1.4 < z < 3.5 
satisfy either sBzK or pBzK criteria. This can be seen 
as either a benefit or a drawback to these selection cri- 
teria; a benefit because only three bands are needed to 
isolate a large number of high redshift galaxies; a draw- 
back because global properties of BzK-selected galaxies 
will be an average over a wide range in redshift, obscur- 
ing the evolving nature of galaxies at an epoch where 
such evolution is expected to be rapid. 

It is also noteworthy that pBzKs have the highest frac- 
tion of galaxies at K < 20 of any of the galaxy popula- 
tions discussed here. The sBzKs in particular have a 
smaller fraction of galaxies this bright at the same red- 
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shifts. This may suggest that pBzKs have a natter lu- 
minosity function in the rest-frame optical than do the 
sBzKs, although a detailed discussion of this point is 
clearly beyond the scope of this paper. A similar con- 
clusion was reached about the luminosity functions of 
z ~ 2.5 red and blue galaxies by Marchesini et al. (2006). 

The photometric redshift distributions of 
BM/BX/LBGs are shown in the bottom panel of 
Figure 12. Excluding the z < 0.4 interlopers, the mean 
photometric redshifts for the three galaxy samples are 
1.4, 2.1, and 2.8, respectively. These values are in good 
agreement with the spectroscopic values of (Steidel et 
al. 2003, 2004), although our if -band selection certainly 
means that our sample will have different properties 
than the bulk of typical 7?.-selected samples (Shapley et 
al. 2004; Adelberger et al. 2005b; Reddy et al. 2005). 
Similarly, the ~ 13% interloper fraction for BM/BX 
galaxies is in excellent agreement with the spectroscopic 
sample of Reddy et al. (2005) (see their Table 3). A 
strength of the BM/BX/LBG selections is that each 
identifies galaxies over a comparatively narrow redshift 
range, so comparisons of the properties of galaxies 
selected by each of these techniques provide meaningful 
constraints on evolution (Adelberger et al. 2005a). 

8.4. Rest-frame optical colors 

Next we investigate the rest-frame optical U-V color of 
galaxies as a function of redshift. The rest-frame colors 
were calculated by interpolating between the observed 
bands using the best-fitting templates as a guide (see 
Rudnick et al. 2003). The results are shown in Figure 13. 
As the J — K color probes the rest-frame optical at z > 2, 
it is no surprise that the DRGs at these redshifts tend 
to be red. The "interlopers" at z < 1.8 have a wider 
range of colors. Compared to the pBzK selection tech- 
nique, the DRG technique does appear to select a larger 
number of red galaxies at z > 2, identifying ~ 75% of all 
galaxies with (U — V) rest > 0.5 (the remaining ~ 25% 
are nearly all at 2 < z < 2.3). For comparison the pBzK 
technique only identifies ~ 30% of these galaxies, but it 
extends to lower redshifts and has very few interlopers. It 
should also be emphasized that the pBzK technique was 
designed to select only passive galaxies, while the DRG 
technique selects both passive and actively star forming 
galaxies (e.g. Labbe et al. 2005: Kriek et al. 2006b); it is 
not currently known if there is significant contamination 
of pBzK samples by star forming galaxies. 

As shown above, the sBzK technique identifies most of 
the galaxies at z > 1.4. Interestingly, there is a signif- 
icant tail of the redshift distribution at 1 < z < 1.4-as 
well as a more sparsely populated tail at lower redshifts 
that seems to be occupied only by bluer galaxies. This 
rcdshift-dependent selection effect may skew any global 
quantities derived from samples of sBzK galaxies. 

The BM/BX/LBG selection criteria identify the bluer 
galaxies. The blue points in the bottom panel of Fig- 
ure 13 meet the < 25.5 limit used for many optical 
surveys, while this magnitude cut has not been imposed 
for the purple points. As expected in a if-selected cata- 
log, galaxies that are fainter in JZ tend to have redder col- 
ors. It has been noted previously that many high redshift 
galaxies do not meet the BM/BX/LBG criteria because 
they are too red-not simply because they are too faint 
in observer's optical (van Dokkum et al. 2004; Daddi et 



al. 2004; van Dokkum et al. 2006). Interestingly, while 
~ 50% of the galaxies at z ~ 2 meet one of the optical se- 
lection criteria, this number increases to ~ 80% at z ~ 3 
(see also the bottom panel of Fig. 12). This suggests that 
a higher fraction of z <~ 3 galaxies have blue continuua 
and prominent Lyman breaks, indicative of significant 
unobscured star formation. This apparent evolution is 
not caused by the Malmquist bias: applying a uniform 
cut in absolute U-band magnitude-the cut is chosen to 
ensures high completeness at z ~ 3-the fraction of galax- 
ies that are selected by any of the BM/BX/LBG criteria 
is ~ 35% and ~ 80% at z ~ 2 and z ~ 3, respectively. 
An important caveat to this result is that z ~ 2 galaxies 
are selected using the BM and BX criteria, while z ~ 3 
galaxies are selected using the LBG criteria. It may thus 
be possible that this effect is an artifact of differences in 
the color selection criteria themselves, rather than evi- 
dence of evolution. However Brammer & van Dokkum 
(2006) study the evolution of optically-red galaxies from 
z ~ 3.7 to z ~ 2.4, and also find evidence for substantial 
evolution in the rest-frame ultra-violet slopes in the same 
sense as described here. 

9. SUMMARY 

We have presented the deep NIR imaging of the 
MUSYC survey. This consists of four 10' x 10' fields, im- 
aged to J - 22.5, H ~ 21.5, and K ~ 21. We combined 
these data with MUSYC optical imaging to produce pub- 
lic if -selected catalogs with uniform UBVRIzJHK pho- 
tometry. 

Many recent surveys rely on a few observed bands to 
isolate large samples of high redshift galaxies. We use 
the high-quality, multi-band photometry from MUSYC 
to investigate some of the properties of galaxies that are 
selected using common selection criteria, including the 
distant red galaxies (DRGs), star- forming BzK galaxies 
(sBzKs), passive BzK galaxies (pBzKs), as well as the 
Lyman break galaxies (LBGs) and the similar BM/BX 
galaxies at somewhat lower redshifts. 

DRGs have a wide distribution of photometric red- 
shifts. The DRGs at z > 1.8 have red rest-frame op- 
tical colors, while the interlopers at lower redshift have 
a range of optical colors. The interlopers account for 
~ 15% of the DRGs at K < 21, but this number increases 
to - 50% at K < 19. In comparison, the BM/BX/LBGs 
each have comparatively narrow redshift windows. If an 
IZab < 25.5 magnitude limit is applied, which is a typ- 
ical limit for ground-based surveys of optically-selected 
galaxies, then the BM/BX/LBGs tend to have blue rest- 
frame optical colors. If this limit is not applied then these 
selection criteria also identify many galaxies with redder 
colors, including some that are DRGs. Interestingly, the 
combined BM/BX/LBG selection criteria identify a sig- 
nificantly higher fraction of the if -selected galaxies at 
z ~ 3 than at z ~ 2; this may suggest an evolution in 
the rest-frame UV properties of red galaxies, as has al- 
ready been found by Brammer & van Dokkum (2006). 
Current samples of DRGs and BM/BX/LBGs are largely 
orthogonal, but the overlap is expected to increase for op- 
tical surveys that reach limits significantly deeper than 
^ab = 25.5. Together, the DRG and BM/BX/LBG cri- 
teria select ~ 90% of the galaxies with 2 < z < 3.5 in 
our sample. 

The sBzK and pBzK criteria select galaxies over a very 
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Fig. 13. — Rest-frame U-V color vs. redshift. The filled circles mark the different galaxy populations discussed in this paper, and 
the small black points mark the remaining MUSYC galaxies. The open purple circles in the bottom panel are for galaxies that meet the 
BM/BX/LBG color criteria, but are fainter than the typical limit of TZab = 25.5 



wide range in redshift. While the sBzK criterion tends to 
select only galaxies with bluer rest-frame optical colors at 
1 < z < 1.4, it also selects many of the reddest galaxies 
at higher redshifts. Because the sBzK criterion selects 
many different types of galaxies over a large range in 
redshift, it may be thought of as a coarse but effective 
photometric redshift technique. Together, the sBzK and 
pBzK criteria detect ~ 93% of the galaxies at 1.4 < 
z < 3.5 to our limit of K = 21. A larger fraction of 
the pBzKs than sBzKs are bright in K, which may be 



indicative of a difference in luminosity functions. Finally, 
the usefulness of the BzK criteria may be limited for the 
reddest galaxies, which may be undetected in even the 
deepest optical images. 

A thorough understanding of the range of proper- 
ties spanned by galaxies selected according to specific 
color criteria requires multiwavelength observations and 
a large sample of spectroscopic redshifts. The use of 
photometric redshifts is a significant limitation of the 
analysis presented in this paper. Unfortunately, a large 
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TABLE 1 
Deep NIR MUSYC Fields 



Field 




RA (J2000.0) 




Dec (J2000.0) 




E(B — V) 


Filter 


— — 

Exposure time (hr) 


HDFS1 


22:33:12.5 


-60:36:40 


0.027 


J 


18.0 










H 


8.3 










K' 


7.0 


HDFS2 


22:31:56.6 


-60:36:46 


0.026 


J 


10.3 










H 


5.0 










K s 


8.5 


1030 


10:30:30.4 


+5:25:00 


0.024 


J 


11.3 










H 


9.9 










K' 


10.4 


1255 


12:55:20.6 


+1:07:49 


0.015 


J 


10.7 










H 


6.9 










K s 


11.3 



number of spectroscopic redshifts for an unbiased sample 
of the galaxies detected by current deep surveys, includ- 
ing those with very red colors, is exceedingly difficult to 
obtain. Progress in the study of galaxy properties as a 
function of redshift will continue to rely heavily on color 
selection techniques, both to define galaxy samples and 
to identify targets for spectroscopic follow-up. An aware- 
ness of the limitations of such techniques is advisable. 
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TABLE 2 

AB CONVERSIONS 



Field filter A eff [A] AB conversion 



HDFS1, 1030 J 12461 0.93 

H 16306 1.38 

K' 21337 1.86 

HDFS2, 1255 J 12470 0.94 

H 16366 1.40 

K s 21537 1.88 



a Defined such that = niYega + conversion 



TABLE 3 

Characteristics of final images 



Field 


Filter 


Zeropoint (Vega) 


FWHM 


5<r depth 


HDFS1 


J 


22.37 ±0.02 


0.96 


22.9 




H 


22.495 ± .015 


0.96 


21.8 




K' 


22.274 ± .02 


0.96 


21.1 


HDFS2 


J 


22.127 ± .025 


1.05 


22.5 




H 


22.522 ± .023 


1.03 


21.4 




K 3 


22.094 ± .022 


1.00 


20.8 


1030 


J 


22.356 ± .02 


0.97 


22.5 




H 


22.545 ± .015 


0.93 


21.8 




K> 


22.23 ± .03 


0.90 


21.3 


1255 


J 


21.959 ± .012 


1.00 


22.6 




H 


22.463 ± .017 


0.92 


21.6 




K s 


22.115 ± .020 


0.93 


21.0 



TABLE 4 

Point-source completeness limits for K-selected 
catalogs 

Field Masking sources Entire image 

90% limit 50% limit 90% limit 5(r%~limit 

HDFS1 21.18 21.50 20.78 21.44 

HDFS2 20.79 21.10 20.55 21.04 

1030 21.37 21.68 21.15 21.59 

1255 20.98 21.30 20.93 21.27 



